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REALIZATION  OF  TARGET  SATELLITE  INTERCEPTION 
WITH  VELOCITY-GAIN  GUIDANCE 

Jia  Peiran  and  Tang  Guojian 

Department  of  Automatic  Control, 

National  University  of  Defense  Technolocpy 


ABSTRACT: 


Xhis  paper  is  about  the  study  of  realization  of  the  target  satellite  interception 
method  with  velocity  gain  guidance.  The  characteristics  are  the  following  j  based  on  the 
given  conditions  of  the  interceptor’s  initial  motion  parameters  and  the  target’s  trajecto¬ 
ry  parameters,  firstly  the  predictive  point  of  impact  intercepting  the  target  under  the  in¬ 
fluence  of  J2  item  of  the  earth  oblatness  is  determined ,  then  the  interceptor’s  velocity  re¬ 
quired  during  the  collision  time  using  Lambert  method  is  decided ,  at  last ,  the  target  in¬ 
terception  with  velocity  gain  guidance  method  is  realized.  A  transform  condition  witjj 
higher  precision  to  the  terminal  guidance  by  using  this  method  is  guaranteed. 

Key  words  interception,  trajectory,  velocity  gain  guidance 


TJ861 

The  article  attempts  to  use  the  velocity-gain  guidance 
method  in  target  satellite  interception  by  an  interceptor. 
Several  problems  are  emphasized  in  the  discussion  of  the  earth 
oblateness  J2  affecting  the  determination  of  the  predicted 
interception  point,  and  the  realization  of  the  guidance  method. 
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In  the  launch  inertial  coordinate  system  determined  by  the 
firing  angle  Aq  and  the  geographic  latitude  Bq,  assume  that  point 
D  in  space  for  the  interceptor  at  time  tj  has  the  following 
motion  parameters ; 

=  [-2:^1)  y^h')  (1) 

VDy(^l)  VD.(il)Y  (2) 

It  is  required  to  locate  the  space  position  point  M  at  which  the 
target  will  be  intercepted  at  tj  +  T  calculated  for  an  orbit 
after  time  T: 

?v(^i  +  T)  —  +  T)  yii/Oi  +  T)  +  T^Y 

correspondingly,  at  times  t^  and  t^  +  T,  the  geocentric  points  of 
the  interceptor  position  at  times  t^  and  ti  +  T  are,  respectively 

?d(^i)  =  +  ^0^  yoW  +  Ro,  ZdUi)  +  R^Y  (4) 

+  T)  =  [xyCi;  +  T)  +  Roi  ytiOi  +  T)  +  R„  Ziti(t\  +  T)  + 

(5) 

In  the  equations,  ^  y  is  the  geocentric-radius  vector 

of  the  firing  point. 

I .  Determination  of  the  Predicted  Interception  Point 

According  to  elliptic  theory,  an  elliptical  trajectory  can 
be  determined  for  two  points  in  space  constrained  in  the  given 
flight  time.  In  addition,  the  required  velocity  at  the  beginning 
point  can  be  determined.  Actually,  the  interceptor  does  not  fly 
in  an  inverse-square  law  gravity  field  law.  After  the  given 
flight  time,  the  interceptor  is  unable  to  collide  with  the 
target,  thus  affecting  the  transformation  condition  for 
transforation  at  the  terminal -guidance  phase.  To  reduce  the 
error  due  to  the  earth's  oblateness,  let  us  adopt  the  method  of 
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the  predicted  interception  point.  After  revising  with  respect  to 
oblateness,  the  interception  point  is  called  the  predicted 
interception  point.  Considering  that  in  the  gravity  field  of  the 
J2  term,  if  the  interceptor  has  the  required  velocity  determined 
by  the  predicted  interception  point,  then  after  flying  the 
interceptor  just  arrives  at  the  required  interception  point. 

First,  the  effect  of  J2  is  not  considered.  According  to 
■ro^rn  and  the  flight  time  T  to  determine  an  elliptical 
trajectory  and  the  instantaneous  required  velocity  S«'. 

(including  velocity  Vr,  velocity  dip  angle  Or,  and  azimuth  angle 
Ar  at  the  elliptic  trajectory  plane).  Moreover,  the 
corresponding  geocentric  latitudes  ro''  and  rn  correspond  to 
points  (|)d  and  <|)m/  as  well  as  the  radial  difference  and 
geocentric  angle  .  Then,  by  using  the  perturbation  theory  [1] 
employed  in  standard  trajectories  to  determine  the  effect  of  the 
J2  term  on  the  interception  point,  and  the  effect  on  flight  time. 
Then  we  can  derive  the  equal-altitude  firing  range  angle  , 

the  equal-altitude  lateral-deviation  angle  Cm^  the  time 

deviation  AA  arrived  at  the  ground  surface,  with  respect  to  the 
distance  between  the  interception  point  and  the  ground  track 
point  that  the  interceptor  arrives  as  determined  by  the 
interception  point  Aj.T.  Convert  Aj-Bm  and  Cm  to  the  linear 
deviation  quantity  of  the  point  of  the  right-handed  inertial 

coordinate  system  of  the  X-axis,  for  the  firing  direction  at  the 
subsatellite  point  with  respect  to  point  D: 


-Axif 

^yu 

At  the  subsatellite  point  for  D,  the  various  deviation  quantities 
in  the  north-sky-east  coordinate  system 
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(7) 


"cosA* 

0 

—  sinA, 

0 

1 

0 

.sinA;{ 

0 

cosAg 

Convert  into  deviation  quantities  in  the  launch  inertial 
coordinate  system 


■/ij/cosfi,  rl,  /ji/cosffel  rAx; 

/ij/cosfB  rly  /jj/cos^b  Ay, 

rl,  /aa/cos^HbJ  LAxj. 


(8) 


In  the  equations. 


/u  =  —  sinforl^ 

fn  =  —  sinfarly 

fu  =  —  sinfork 

f  1\  —  ^z^Dy 

fz2  =  Oyo.  -  /3»rL 
f»  =  ii>ly  -  Olrly 
VL”'  r^ox  + 

^or  +  yo 
riz  Li?„.  4- Zd- 

01  rcos5(,cosAo 

OPy  =  sin£, 

Plj  L—  cos5„sinAo- 


(9) 


(10) 


(11) 


Finally,  reduce  the  deviation  quantity  due  to  oblateness  from  the 
interception  point  ,  then  the  position  parameter  :ri  of  the 
predicted  interception  point  is 

ri  =  [r^ix  ti«]^  =  [^ox  +  +  yit  ^oz  +  (12) 
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In  the  equation, 


Xv  ~  -As:' 

v.i 

yn  — 

Xu  — 

(13) 


There  is  a  time  deviation  A^T  as  the  interceptor  flies  to 
the  interception  point  in  the  inverse-square  law  gravity  field. 
Therefore,  a  revision  should  also  be  made.  In  other  words,  it  is 
required  that  the  interceptor  fly  to  rm  in  the  inverse-square 
law  gravity  field  within  the  time  T*=T-ArT. 

II.  Determination  of  Required  Velocity  for  Fixed  Interception 
Time 


Now  the  problem,  based  on  the  real-time  position  ri  and  the 
position  ro'  of  the  predicted  interception  point  of  the 
interceptor,  is  to  determine  the  required  velocity  vj  for  the 
interceptor  to  arrive  at  the  predicted  interception  point  in  an 
inverse-square  gravity  field  within  the  time  period  T*.  This 
point  can  be  solved  by  using  the  Lambert  method  [3].  Moreover, 
from  the  geocentric  vector-radii  rc>  and  ri:  ,  we  can  determine 
the  sine  and  cosine  values  of  the  great-arc  azimuth  angle  AR  of 
the  elliptical  orbital  plane  defined  by  ^d'  and  ;ri  ,  from  the 
longitude  difference  /lo«:  .  between  the  two  points  of  geocentric 
latitudes  (t)^  and  (pM. 

sin>i;  =  cos?*  sinl^v/sin/Si 

'  =  (cos?*  -  cosy9isin?fc)/(sin/?*cos?b)  (24) 

./?i  =  arccos[?„  .  r*/(r£,ri)] 


Thus,  it  is  required  that  the  projection  of  velocity  vR  in  the 
launch  inertial  coordinate  system  be  given  as 


'vL' 

7u/cos?b  rL  fijcosfo' 

cos0^  coSi4j; ' 

—  Vr 

/12/cos?®  r%y  f%i/cos<po 

sinSj? 

_cos0;  sinAg  _ 
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In  the  equation,  0-  is  the  dip  angle  of  the  required  velocity 

Vk  . 

III.  Realization  of  Velocity-Gain  Guidance  Method 

Generally,  we  call  the  difference  between  the  instantaneous 
required  velocity  v'r  and  its  real-time  velocity  I'ci  of  the 
interceptor  as  follows ; 

V,  =  Vr  —  vd  (16) 

In  the  equation,  5,  is  the  velocity  gain. 

In  principle,  align  the  booster-acceleration  direction  and 
the  gain-velocity  direction  of  the  interceptor  so  that  after  time 
tau  of  engine  burn,  the  interceptor  measurement  apparatus 
provides  velocity  and  position  parameters  at  a  new  time  instant, 
and  repeat  the  foregoing  steps  to  continue  the  guidance.  We 
should  point  out  that  if  the  time  is  T  as  specified  for  the 
interceptor  from  the  beginning  guidance  point  to  the  interception 
point,  then  n  iterations  will  subtract  n  times  tau  from  T. 

Besides  the  foregoing  fundamental  concept,  the  following 
problems  can  be  considered  in  engineering  practice. 

3.1.  Estimation  of  required  velocity,  and  improvement  of 
velocity-gain  at  engine  shutoff  point 

To  save  fuel  and  to  allow  the  interceptor  have  a  more  stable 
attitude  near  the  engine  shutoff  point,  we  consider  using  the 
required  velocity  vk.k>  at  the  estimated  engine  shutoff  point  to 
replace  the  instantaneous  required  velocity  ^1.. 

Series  expansion  at  the  time  instant  t^  with  respect  to  the 
velocity-gain  v  and  the  required  velocity  'vk  ,  we  select  the 
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(17) 

(18) 


following,  in  approximate  terms; 


~  ^,(i*)  =  w,(f,)  +  ?,(f,)(fj  —  ;.) 
vlt  =  (fj  =  5;  (f.) 


In  the  equation. 


as  lw,(A)  -  ?,(«._, )]/r 
t'x  (<.)  *=  [Sr  (f.)  —  Sr  (t._,)]/r 
r  =  f.  -  f.-i 


we  also  know  when  tj^ - >t^^,  vg,k - >0,  therefore  from  Eg.  (17), 

we  know 


‘i-ti  =  —  [S,(ti)  •  1  v,(ti)  i* 

S;.i  =  S«*  (is,)  -  -^[Sr*  (<,)  -  vi  (/.-i)][S,(i,)  •  S,(r,)]/!  S^Ct,.) 


(19) 

(20) 


Eg.  (20)  is  the  estimated  vector-equation  for  the  velocity- 
required  at  the  engine  shutoff  point. 

Consideration  is  given  to  the  effect  of  gravity  in  the  time 
interval  between  t^  and  t^,,  and  the  gravity  acceleration  within 
this  time  period  is  considered  as  the  to  ,  then  the  improved 
velocity-gain  of  the  interceptor  can  be  expressed  as: 


=Sr  (;.)  — SD(it)  —  — ?i>(t.)  [S,(r,)  •  i;,(;,)]/|i!,(f,)|*  (21) 


3.2.  Determination  of  interceptor  attitude 
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Assvime  that  the  booster  apparatus  of  the  interceptor  is 
composed  of  four  constant-value  small-booster  engines  vertical  to 
the  interceptor  main  axis  o^Xj,  as  installed  on  the  OiYi  and  the 
OiXi  axis  cruciformly  surrounding  the  mass  center.  When  applying 
thrust-vectoring  fixation  in  the  v,'  direction,  the  attitude- 
control  scheme  of  the  interceptor  is  various.  If  we  base  our 
calculations  on  the  principle  of  minimum  change  in  attitude 
angle,  to  select  the  single  Oiyi  and  the  OjZi  axis  direction  of 
the  engine,  in  the  v/  vector  direction,  to  provide  thrust,  or 
to  allow  the  combined  thrust  of  the  oiy^  and  OjZj  axes  with  v,'  to 
be  coplanar  to  simultaneously  start  two  engines  for  guidance 
control.  A  simpler  scheme  is  to  use  the  positive  thrust  from  the 
OiYi  axis  to  align  the  direction  with  the  -5,  direction,  for 
guidance.  If  it  is  assumed  that  the  interceptor  coordinate 
system  at  the  beginning  guidance  point  is  coincident  with  the 
platform  coordinate  system,  then  the  attitude  control  angle  of 
the  interceptor  is 

f 

'  'P 

y 


=  arctg(L',,/t;',,)  —  sign(x-,^  •  y) ,  —  n-  <  f  <  ;r 

=  0  (22) 
=  arctg(t;^/  “  J  ^  f- 


3.3.  Control  of  engine  shutoff  point 

When  conducting  computations  for  the  engine  shutoff 
condition,  the  interceptor-borne  computer  has  a  time  lag  error. 
Assume  that  the  computation  step  length  is  t,  then  the  engine 
shutoff  time  error  can  be  between  tau  and  2t.  To  reduce  the 
control  error  at  the  engine  shutoff  point,  the  two  following 
measures  can  be  considered  to  upgrade  guidance  accuracy. 

(1)  A  linear  prediction  is  conducted  on  the  engine  shutoff 
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time  in  order  to  determine  the  time  of  simplified  computation 
using  the  guidance  equation. 

Vf-  at  time  t  near  the  engine  shutoff  point  can  be  considered 
as  a  linear  function  of  time.  Then,  Eq.  (19)  can  be  used  to 
determine  the  time  interval  (t^^-t)  from  t  to  to  to  tj^,,  the  time 
r  of  engine  shutoff  point.  Allowance  is  given  to  the  fact  that 
before  the  use  of  the  small-step  length,  the  time  length  error  of 
the  interceptor-borne  computations  is  r'  .,  the  original 
computation  step  length.  Therefore,  selecting  the  time  for 
simplified  computations  of  the  guidance  equation  is  the  first  to 
satisfy  the  time  t 

^  <  2r  +  2r  (23) 

(2)  Simplified  equation  of  guidance  equation  when 
calculating  small-step  length  near  engine  shutoff  point 

After  determining  the  time  t  satisfying  Eq.  (23),  actually 
the  interceptor-borne  engine  has  operated  for  a  step  t.  By 
marking  tx=t+T,  thus  beginning  at  tT,  the  small-step  length 
computation  commences.  Generally  speaking,  only  when  the 
computer  run  time  does  not  exceed  t+2t,  which  is  a  very  short 
time,  therefore,  the  general  guidance  equation  can  be  simplified. 
If  the  time  parameter  T  is  used  as  the  basis,  the  real-time 
velocity  is  5c'.  ,  the  required  velocity  is  vi.  ,  then  the 
extrapolated  velocity  gain  5,'  of  the  interceptor  can  be 
obtained . 

VdUt  +  )  =  ?c(<)  -f  —Ivod')  —  v(t  —  r)](r  +  nr")  (24) 

vi  (h  +  riT  ')  =  vi  (t)  4-  (O  —  viCt  —  r)](r  +  nr)  (25) 
v,(tT  +  «r')  =  ?;  (fr  +  rtf')  —  VoOt  +  nr')  (26) 
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In  the  equations,  n=l,2,...  is  the  number  of  iterative 
computations  after  using  the  small-step  length.  Thus,  Eq.  (19) 
can  be  used  to  calculate  the  value  (tk-tT-nx').  Take  note  that 
guidance  of  the  interceptor  is  not  necessary  within  a  very  short 
time  near  the  engine  shutoff  point.  Now,  the  attitude  is  the 
attitude  of  the  final  large-step  length.  The  reason  is  that  the 
tiny  variation  in  the  interceptor  attitude  within  a  very  short 
time  does  not  have  much  effect  on  the  motion  of  the  mass  center. 
Besides,  the  slowly- varying  gu  is  considered  as  a  constant.  By 
considering  the  accuracy  of  the  linear  prediction  and  the 
feasibility  of  prediction  when  calculating  the  time  delay  and  the 
engine  shutoff  time,  the  first-appearing 


~  ir  —  nr'  ^  2*^ 


(27) 


is  considered  as  the  determination  condition  of  issuing  the 
linear  prediction.  Then,  when  the  above  equation  is  established, 
the  computation  formula  for  the  engine  shutoff  time  is 


h  —  tr  nr' 


(28) 


IV.  Numerical  Simulation  Results  and  Discussion 

4.1.  Simulation  results 

Assiime  that  the  constant  thrust  acceleration  generated  by 
each  engine  is  O.Sgo,  when  the  foregoing  scheme  is  adopted,  the 
simulation  results  of  a  typical  interception  trajectory  is  as 
follows;  the  guidance  operating  time  is  29.196s,  at  the  150km 
distance  between  the  interceptor  and  the  target  satellite,  when 
the  terminal  guidance  phase  begins,  the  line-of-sight  rotating 
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rate  for  the  relative  motion  of  the  two  is  3 . 096xl0"®/s ;  this 
entirely  satisfies  the  initial  conditions  for  terminal  guidance. 

4.2.  Analysis  and  discussion 

(1)  By  using  the  predicted  interception  point  after 
modification  with  the  oblateness  J2  term,  the  velocity-gain 
guidance  is  conducted,  along  with  measures  of  simplified 
computation  at  the  predicted  engine  shutoff  time  and  near  the 
engine  shutoff  point,  along  with  the  use  of  small-step  length 
computation,  the  error  is  very  small  when  use  is  made  of  the 
method  of  guidance  based  on  simulation  results .  When  there  is  no 
external  interference,  in  the  guidance  process  the  engine  is 
started  up  only  once.  It  is  apparently  that  the  use  of  the 
velocity-gain  guidance  method  to  accomplish  the  guidance  task  is 
feasible . 

(2)  In  engineering  practice,  if  problems  of  computer  memory 
and  operating  speed  are  encountered  when  using  this  guidance 
method,  consideration  is  given  not  to  use  the  method  of  the  post- 
modified  predicted  intersection  point,  but  to  use  the  required 
velocity  vR  for  guidance  in  the  inverse-square  law  gravity 
field  obtained  from  the  given  interception  point.  In  other 
words,  multiple  engine  firings  are  required  to  continuously 
update  the  error  in  guidance  when  not  considering  the  effect  of 
oblateness,  to  improve  the  guidance  accuracy. 


The  article  was  received  for  publication  on  May  2,  1991. 
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